Ginzburg-Landau- Wilson Aspects of De-confined Criticality 
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Monte Carlo study of the deconfined critical action phase diagram reveals a region where spinon 
deconfinement occurs through a weak first-order phase transition, in agreement with Ginzburg- 
Landau theory. Wilson renormalization argument in combination with the absence of the data 
collapse even in the regime of weak interaction between the spinons casts a serious doubt on the 
possibility of the continuous deconfinement transition. We also argue that if a continuous deconfined 
criticality does exist on the phase diagram, its nature is analogous, in a certain precise sense, to 
that of the XY universality class. 



CO 



I 

C 

O 

o 



(N 
> 
(N 

in 
O 

o 
in 
o 



-a 
c 

o 
o 



X 



PACS numbers: 05.30.-d,75.10.-b,05.50.+q 

Critical properties of systems described by one or sev- 
eral complex fields coupled to one gauge field have a 
long history of studies and apply to numerous problems 
in physics which include normal-superfluid transitions 
in multi-component neutral or charged liquids (see, e.g. 
0), superfluid-valence-bond solid (SF-VBS) transitions 
in lattice models (0,0, H)> Higgs mechanism inparticle 
physics etc. Recently, the authors of Refs. 0,0 ar- 
gued that the SF-VBS transition in a (2 + l)-dimensional 
system is an example of a qualitatively new type of quan- 
tum criticality that does not fit the Ginzburg-Landau- 
Wilson (GLW) paradigm. The discussion of this claim is 
the main focus of our Letter. 

The theory of deconfined critical point (DCP) estab- 
lishes a remarkable microscopic picture of how a generic 
continuous SF-VBS transition may happen. A generic II- 
order SF-VBS transition can not be derived from naive 
GLW expansion in powers of order parameters starting 
from the quantum disordered groundstate, which has no 
broken symmetries, simply because such a state is un- 
likely to be present in most models 0, @ ■ It may hap- 
pen, however, that restrictions imposed on parameters of 
the GLW action are automatic and originate from hidden 
(or emerging) symmetries. In the DCP action the cru- 
cial symmetry of this kind is that between the spinons (or 
vortices) in the VBS state. Another alternative is that 
the DCP action in fact predicts that the SF-VBS phase 
transition is of the I-order. 

In what follows we discuss results of the Worm algo- 
rithm Monte Carlo simulations performed for the DCP 
action, construct its phase diagram, and show that the 
SF-VBS transition is I-order at least for sufficiently 
strong coupling between the spinon fields. This explains 
the outcome of recent numerical simulations which re- 
ort either difficulties with finite-size scaling of the data 
or weak I-order transitions |g. It is far more diffi- 
cult to present a direct evidence of the I-order transi- 
tion in the weak coupling limit. However, problems with 
the data collapse in finite-size systems hint at this pos- 
sibility. The global I-order scenario would not be a sur- 
prise in view of a simple renormalization argument that 



in the DCP action the weak coupling regime maps onto 
the intermediate-coupling one, where we find a circum- 
stantial evidence that the Il-order scenario fails. 

Many features of the DCP action are remarkably sim- 
ilar, both qualitatively and quantitatively, to those of a 
more conventional two-component XY-model (2XY) 
Moreover, we argue that if the line of Il-order DCP's does 
exist, then its critical properties are most close to those 
of 2XY at the U(l)xU(l) critical point, in the follow- 
ing sense. There exists a self-dual model with marginal 
long-range interactions oc 1/r 2 , which continuously in- 
terpolates between DCP and U(l)xU(l) criticality by 
varying amplitudes of oc l/r 2 terms. 

There are many equivalent formulations of the action 
describing coupled gauge and multi-component complex 
fields 0, 0, @ . Here we employ the integer-current lat- 
tice representation which can be viewed either as a high- 
temperature expansion for the XY model in three dimen- 
sions, or as a path-integral (world-line) representation 
of the interacting quantum system in discrete imaginary 
time in (d + 1) =3 dimensions. The XY action reads 



U T — r'jr ' Jr 
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where summation over repeated lattice sites r is assumed; 
U T - r > = U5 Tt r'', and j r = (j r )^ with // = x,y,T are inte- 
ger, zero-divergence, Vj = 0, currents defined on bonds of 
the simple cubic space-time lattice with periodic bound- 
ary conditions. The configuration space of j-currents is 
that of closed oriented loops. In terms of particle world 
lines, currents in the time direction represent occupation 
number fluctuations away from the commensurate filling, 
and currents in the spatial directions represent hopping 
events. The state of S with small current loops is normal 
("high-temperature"); the corresponding particle order- 
parameter, ip, is zero in this phase. When particle world 
lines proliferate and grow macroscopically large, the sys- 
tem enters into the superfluid state with (ip) ^ 0. 

Generalization to the symmetric multi-component case 
with short-range interactions is straightforward (in what 
follows we will concentrate on the two-component case 
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FIG. 1: Phase diagrams of the long-range (left panel), Eq. Q, 
and short-range (right panel), Eq. JSJ, actions. Bold solid 
lines indicate regions where I-order phase transitions were un- 
ambiguously verified. Error bars are shown but are typically 
much smaller than symbol sizes. 
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FIG. 2: Normalized probability density distributions P(S) 
of the DCP action for g = 1.5 and different system sizes; 
U(L=8)=0.94375, U(L=12)=0.92345, U(L=16)=0.91667, 
U(L=22)=0.91274. 



which describes coupled XY models) 

S { 2 S) = Ujj r + U&-V (j lr + j 2r ) 2 . (2) 

If particle/XY-spin fields in the original model are cou- 
pled through a gauge field, then after the gauge field is 
integrated out one arrives at the DCP action similar to 
Eq. which now contains a long-range part 0, E| 

sf = Uf lr + Uf 2r +gQ M , (j lr + j 2r ) •(jir' + j 2r <) . (3) 

The lattice Fourier transform of the interaction potential 
Qr-r' is given by Q q = 1/ ^2 sin 2 (q A1 /2), which implies 
the asymptotic form Q(r — > oo) ~ 1/r. In the discussion 
of the SF-VBS transition, jx and j 2 currents in the DCP 
action represent world lines of spinons which are VBS 
vortices carrying fractional particle charge ±1/2 |2|,|3j. 

Phase diagrams for the long- and short-range actions 
are shown in Fig. ^ Different phases are identified in 
terms of loop sizes for j\ and j 2 currents as follows: the 
Mott insulator (MI) and VBS states are characterized 
by small j\- and j 2 -loops; in the VBS supersolid (SFS) 
and paired superfluid (PSF) states only (jx =p j 2 )-loops 
grow macroscopically large while single component loops 
remain small; finally, in the SF state there are macroscop- 
ically large j\- and j 2 -loops. The data for the short-range 
model are reproduced from Ref. S|- A similar shape of 
the phase diagram for the DCP action in the angle-gauge 
field representation was obtained in Ref. ^(j i though, no 
I-order transition was identified [Tl| . 

It is hard not to notice a remarkable quantitative sim- 
ilarity between the phase diagrams. In particular, the 
DCP action has VBS to SFS and SFS-SF transitions 
when coupling between the spinons is strong and they 
form tightly bound particle states. One immediate pre- 
diction based on properties of the short-range model 



would be that close to the SF-SFS-VBS bicritical point 
the SF-VBS transition is I-order in nature. Unfortu- 

(s) 

nately, the I-order transition in S 2 is very weak and can 
be clearly seen only at system sizes L > 24 by calculating 
the probability density for the system action, P(S), and 
observing its double-peak structure. Simulating large L 
for the long-range action is far more difficult. We were 
able to collect reliable statistics for S 2 on ly f° r L < 22. 
Still, for g = 1.5 we resolved (Fig. [2J the development 
of the double-peak structure in P(S) starting from an 
anomalously flat maximum at L = 8. 

For small g one expects the I-order transition to be- 
come weaker (the g — point describes two indepen- 
dent XY models) and nearly impossible to study us- 
ing P(S) functions. However, there are two important 
circumstances, (i) At g — > 0, the initial part of the 
renormalization flow with L — > oo simply leads to in- 
creasing g(L) oc L 01 [ see aLS0 Eq. ©], thus mapping 
the weak-coupling regime onto the intermediate-coupling 
one. [Indeed, 1/r interaction is a relevant operator re- 
sponsible for the spinon confinement.] (ii) For the su- 
perfluid stiffness, n s , a scale invariant Il-order critical- 
ity implies n s L — > const as L — > oo. These facts has 
to be combined with the numerical observation that for 
g > 0.1 the n s L curves do not intersect at the same U 
even approximately (see Fig. 0) . This strongly suggests 
the I-order transition scenario for all g (we exclude a 
non-scale-invariant Il-order scenario). 

Critical self-duality of XY models. Speaking strictly, 
one cannot rule out that a Il-order line actually starts 
below g = 0.1, or, that for some reason the data collapse 
for g = 0.5 sets in only at L > 20. A possibility of Il- 
order phase transitions in generic multi-component XY 
models with 1/r current-current interactions can hardly 
be questioned since the duality transformation for 
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the short-ranged action S 2 produces a family of such 
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FIG. 3: n s (U)L curves for g = 0.5 and system sizes L 
6,8,12,16,22. 



models. To be unambiguous in our further analysis, we 
adopt the following two definitions. 

Def. 1: We call a model self-dual if there exists a duality 
transformation which at the critical point preserves the 
form of the Hamiltonian, but not necessarily the values 
of the critical parameters. 

Def. 2: A model A is critically self-dual if there exists a 
self-dual model B such that the critical properties of A 
and B are the same. 

XY action can be identically transformed into the dual 
(inverted- XY) action which has the same functional 
form as Eq. Q for integer, zero-divergence currents 1 
coupled by the long-range dual interaction potential, 
[7r5oo ~ 1/r. In Fourier space, the dual potential is 
given by Uq = (n 2 /4)Q q /£7 . The notion of bond cur- 
rents also changes from particle to vortex world lines. 
Now, the proliferation of vortex 1-lines signals the onset 
of the transition to the normal state. Formally, long- 
range interactions in the dual action make the XY model 
not self-dual in conventional sense, though the configura- 
tion space is preserved. Several papers suggested that 1/r 
interactions result in the inverted- XY criticality qualita- 
tively different from the conventional XY point [lj, [l5| . 

A single-component short-range XY model is criti- 
cally self-dual. Below we prove this statement numer- 
ically. Meanwhile, it is easy to show analytically that 
XY criticality is arbitrarily close to the critically of self- 
dual models. First, we introduce XY models with 1/r 2 - 
potential, which are self-dual. Indeed, the potential 
with the Fourier transform i7 a _»o 



A(Tr/2)y/Q(g) 



under duality transformation becomes U t 



(d) 
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A~ 1 (n/2) y /Q(q) - B/A 2 + .... Note that the amplitude 
in front of the 1/r 2 term changes. Second, we observe 
that the 1 /r 2 interaction is marginal (this will be explic- 
itly seen below) and thus by adding a l/r 2 -term with 
infinitely small amplitude to XY model JTJ we do not 
change its critical behavior and exponents, but formally 
make it self-dual. It is of crucial importance that XY crit- 
icality is mapped onto that of a model with marginal in- 



teraction; this is required to reconcile critical self-duality 
with different critical exponents and universal numbers 
of the XY and inverted XY transitions. These differ- 
ences naturally follow from different amplitudes of the 
dual l/r 2 -terms. The XY and inverted XY critical points 
turn out to be just two points in the continuum of A/r 2 - 
criticalities, corresponding to two special values, A^ and 
A + . By changing A one may continuously interpolate be- 
tween A- and A + , gradually transforming XY criticality 
into inverted XY, and vice versa. 

To illustrate the above-mentioned point and, in par- 
ticular, to extend the proof to the case of finite A, we 
computed the critical exponent Dh (A) for the interaction 
potential t/ q = A(ir/2)^Q(q)+B. The Hausdorff dimen- 
sion Dh gives the the average line length of critical loops 
as a function of system size, (I) cx L Dh . Under dual- 
ity transformation. Dh (A) transforms into D^ (A) . As 
expected, among a continuum of self-dual critical points 
B C (A) there is one, A « 3.2, B C (A = 3.2) = -2.8457(10) 
which reproduces Hausdorff dimensions of the inverted- 
XY and XY systems respectively within the statistical 
error bars jl6| . 

To understand the origin of critical self-duality, and in 
particular, the physical meaning of the 1/r 2 interaction, 
we start with constructing proper coarse-grained current 
variables similar to winding numbers in finite-size sys- 
tems and subject to the real space-time renormalization 
group (RG) transformations. Crucial for our derivation 
will be the fact that at criticality there are ~ 0(1) loops 
of size ~ r in the volume r 3 |17| . which is nothing but 
the requirement of scale invariance. 

Imagine a cube of linear size u centered at point r, 
where u 3> 1 is odd integer. Its facets cut bonds of the 
original lattice a distance u/2 away from the center. In- 
teger, zero-divergence currents J r are defined as current 
fluxes through the facets of the cube (we denote the cube 
by C r and its facets as S r:fl ) 



(Jr) M = 



(4) 



r>£S r 



with an additional restriction that only j-loops of size 
larger than u/2 contribute to the sum Q. The last con- 
dition is necessary to suppress noise due to small loops 
winding around the cube edges. The new variables are 
scale-invariant by construction, i.e. (J 2 ) ~ 0(1). This 
fact alone allows us to study the RG flow of the long- 
range effective interaction potential. Let the original 
model be long-range with the bare potential decaying as 
g/r a , a<2. After coarse- graining the interaction energy 
between two cubes C ri and C T2 separated by a distance 
R = |ri — r 2 | » u can be estimated as 



e{u)R c 



E E j> 

xec ri yec r2 



Jy 



9 u 



e{u)R< 



- Jri ' Jr 2 7 



(5) 



where e(u) is the "dielectric" constant at length-scale u. 
The last relation follows from the continuity of world 
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lines and the fact that J-currents are facet, not bulk, 
sums, so that X)xec J x ~ u ^ri- We now rescale all 
distances by a factor of u and observe that the effective 
interaction potential transforms as gu 2 ~ a / e(u)r a . This 
means that (i) long-range interactions with a < 2 are 
screened at the fixed point, e(u) oc u 2 ~ a — > 0, and (ii) 
the term J r • J r '/( r — r ') 2 i s a marginal operator. For 
q = 1 and small initial values of g, the RG flow starts 
with epsilon = 1 and leads to g(u) ~ u. This result 
was used to relate the weak- and intermidiate-coupling 
regimes in the DCP action. This consideration does not 
depend on the number of components, but it works only 
as long as g{u) remains small. 

The above analysis is supported by microscopic mecha- 
nisms leading to the renormalized l/r 2 -law between mi- 
croscopic currents. In the SF phase, the 1/r potential 
originates from the kinetic energy of the superflow around 
vortices. Close to the critical point, the superflow around 
a vortex line is screened on large distances by small vor- 
tex loops, and the kinetic energy integral is transformed 
into (for simplicity we consider here a straight vortex 
line) J n s (r)dr/r — > J dr/r 2 with scale-dependent super- 
fluid density n s oc 1/r. This formula implies that at the 
critical point the interaction energy between the vortex 
line elements is proportional to 1/r 2 . The emergence of 
the renormalized potential ~ 1 jr 2 while approaching the 
critical point from the normal phase, can be considered 
as a result of virtual exchange of long-wave sound exci- 
tations. Elementary calculations in the hydrodynamical 
approximation for d — 2 superfluid show that such po- 



tential scales as ~ 1 jr 2 . 

One may also take an alternative, mathematically rig- 
orous, point of view and consider an action describing 
a single j-current line (of arbitrary fractal structure) at 
the critical point when other degrees of freedom in the 
system are integrated out 



S. 



(loop) 



(^crit)r-r' + A ^r-r 



Jr ' Jr' 



(6) 



where f7 cr it = A/r 2 + short-range terms. It is written in 
terms of original j- or /-currents and not coarse-grained 
variables. In particular, the statistics of the discon- 
nected loop, formally identical to a polymer of a variable 
length, must reproduce geometrical exponents at the crit- 
ical point. It is this critical-loop action that unifies short- 
and long-range XY models, putting them in a wider con- 
text of self-dual 1/r 2 models. 

In conclusion, we present numerical and RG agruments 
strongly suggesting that the DCP action is the theory of 
weak I-order SF-S transitions. In the "zoo" of self-dual 
1/r 2 criticalities interpolating between short-range and 
long-range 1 jr models, the number of components does 
not seem to play any qualitative role unless the transition 
turns into the I-order one. 
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